Grain Boundary Scars and Spherical Crystallography 
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We describe experimental investigations of the structure of two-dimensional spher- 
ical crystals. The crystals, formed by beads self-assembled on water droplets in oil, 
serve as model systems for exploring very general theories about the minimum energy 
configurations of particles with arbitrary repulsive interactions on curved surfaces. 
Above a critical system size we find that crystals develop distinctive high-angle grain 
boundaries, or scars, not found in planar crystals. The number of excess defects in a 
scar is shown to grow linearly with the dimensionless system size. The observed slope 
is expected to be universal, independent of the microscopic potential. 



Spherical particles on a flat surface pack most effi- 
ciently in a simple lattice of triangles, similar to billiard 
balls at the start of a game. Such six-fold coordinated 
triangular lattices Q cannot, however, be wrapped on 
the curved surface of a sphere; instead, there must be 
extra defects in coordination number. Soccer balls and 
Ceo fullerenes 0, 3J provide familiar realizations of this 
fact - they have 12 pentagonal panels and 20 hexago- 
nal panels. The necessary packing defects can be char- 
acterized by their topological or disclination charge, q, 
which is the departure of their coordination number c 
from the preferred flat space value of 6 (q = 6 — c); a 
classic theorem of Euler 0, 0] shows that the total discli- 
nation charge of any triangulation of the sphere must be 
12 [6(. A total disclination charge of 12 can be achieved 
in many ways, however, which makes the determination 
of the minimum energy configuration of repulsive par- 
ticles, essential for crystallography on a sphere, an ex- 
tremely difficult problem. This was recognized nearly 
100 years ago by J.J. Thomson 7], who attempted, un- 
successfully, to explain the periodic table in terms of rigid 
electron shells. Similar problems recur in fields as diverse 
as multi-electron bubbles in superfluid helium ||, virus 
morphology 11(1 Hlj . protein s-layers and cod- 

ing theory [14i llq. Indeed, both the classic Thomson 
problem, which deals with particles interacting through 
the Coulomb potential, and its generalization to other in- 
teraction potentials remain largely unsolved after almost 
100 years [3E1E3- 

The spatial curvature encountered in curved geome- 
tries adds a fundamentally new ingredient to crystallog- 
raphy, not found in the study of order in spatially flat sys- 
tems. To date, however, studies of the Thomson and re- 
lated problems have been limited to theory and computer 
simulation. As the number of particles on the sphere 
grows, isolated charge 1 defects are predicted to induce 



too much strain; this can be relieved by introducing ad- 
ditional dislocations, consisting of pairs of tightly bound 
5-7 defects [l9| which still satisfy Euler's theorem since 
their net disclination charge is zero. Dislocations, which 
are point-like topological defects in two dimensions, dis- 
rupt the translational order of the crystalline phase but 
are less disruptive of orientational order . While they 
play an essential role in crystallography on a spherical 
surface, the configuration and orientation of these ex- 
cess defects remains undetermined, and can only be fully 
understood through a combination of theory and exper- 
iment. Experimental realizations that probe the subtle 
structures have, however, been sorely lacking. 

We present an experimental realization of the gener- 
alized Thomson problem that allows us to explore the 
lowest energy configuration of the dense packing of re- 
pulsive particles on a spherical surface. We create two 
dimensional packings of colloidal particles on the sur- 
face of spherical water droplets, and view the structures 
with optical microscopy. Above a critical system size, the 
thermally equilibrated colloidal crystals display distinc- 
tive high-angle grain boundaries, which we label "scars" . 
These grain boundaries are found to end within the crys- 
tal, which is not observed to occur on flat surfaces be- 
cause the energy penalty is too high. 

Our experimental system is based on the self-assembly 
of one micron diameter cross-linked polystyrene beads 
adsorbed on the surface of spherical water droplets (of 
radius R), themselves suspended in a density-matched 
toluene/chlorobenzene mixture |20j . The particles are 
imaged with phase contrast using an inverted microscope. 
The curvature of the spherical water droplet limits the 
imaged surface area to between 5% and 20% of the full 
surface area of the sphere, depending on the size of the 
droplet. After determining the center of mass of each 
bead, the lattice geometry is analyzed by Delaunay tri- 
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angulation algorithms |2lJ appropriate to spherical sur- 
faces. 

We analyze the lattice configurations of a collection 
of 40 droplets. A typical small spherical droplet with 
system size, R/a = 4.2, where a is the mean particle 
spacing, is shown in Fig. The associated Delaunay 
triangulation is shown in Fig. 03. The only defect is one 
isolated charge +1 disclination. Extrapolation to the en- 
tire surface of the sphere is statistically consistent with 
the required 12 total disclinations. 

Qualitatively different results are observed for larger 
droplet sizes as defect configurations with excess dislo- 
cations appear. Although some of these excess dislo- 
cations are isolated, most occur in the form of distinc- 
tive (5 — 7 — 5 — 7 — • • • — 5) chains, each of net charge 
+1, as shown in Fig. QD. These chains form high-angle 
(30°) grain boundaries, or scars, which terminate freely 
within the crystal. Such a feature is energetically pro- 
hibitive in equilibrium crystals in flat space. Thus, al- 
though grain boundaries are a common feature of 2D 
and 3D crystalline materials, arising from a mismatch 
of crystallographic orientations across a boundary, they 
usually terminate at the boundary of the sample in flat 
space because of the excessive strain energy associated 
with isolated terminal disclinations. Termination within 
the crystal is a feature unique to curved space. Thus, 
our results provide important guidance to determine the 
configuration of excess defects on a sphere. 

To ascertain that the colloidal particles are equili- 
brated, we use particle tracking routines and subsequent 
automated triangulation to measure the mobility of the 
observed screening dislocations. The diffusion of ther- 
mally excited colloid particles on the surface of the wa- 
ter droplets results in local rearrangements of the crystal 
structure. Thermal fluctuations create and destroy dis- 
locations once every few seconds, on average, indicating 
that the defect arrays reach equilibrium much faster than 
the observation time of 10-60 min. The equilibration time 
can also be estimated by the time required for a dislo- 
cation to diffuse across typical defect structures. The 
measured diffusion constant allows us to calculate equili- 
bration times that range from a few seconds to hundreds 
of seconds, depending on the size of the crystal. Since the 
spherical crystals exist for 10 to 60 minutes, the system 
has sufficient time to reach equilibrium. Thus, our obser- 
vations reflect the equilibrium ground state, as opposed 
to a history-dependent nonequilibrium effect, which is 
the case for crystals in flat space. 

To quantify the behavior of the scars, we determine the 
number of excess dislocations per chain for each droplet, 
with the convention that dislocations are counted as be- 
ing part of the same array if they are within three lat- 
tice spacings, and plot the results as a function of R/a 
in Fig. El Scars only appear for droplets with R/a > 5. 
These results provide a critical confirmation of a theoret- 
ical prediction that R/a must exceed a threshold value 




FIG. 1: Light microscope images of particle-coated droplets. 
Two droplets (A) and (C) are shown, together with their 
associated defect structures (B) and (D). Panel (A) shows 
an ~ 13% portion of a small spherical droplet with radius 
R = 12.0 microns and mean particle spacing a = 2.9 mi- 
crons (R/a = 4.2), along with the associated triangulation 
(B). Charge +1( — 1) disclinations are shown in red and yel- 
low respectively. Only one +1 disclination is seen. Panel (C) 
shows a cap of spherical colloidal crystal on a water droplet of 
radius R — 43.9 microns with mean particle spacing a — 3.1 
microns (R/a = 14.3), along with the associated triangulation 
(D). In this case the imaged crystal covers about 17% of the 
surface area of the sphere. The scale bars in (A) and (C) are 
5 microns. 




FIG. 2: Excess dislocations as a function of system size. The 
number of excess dislocations per minimal disclination iV as a 
function of system size R/a, with the linear prediction given 
by theory shown as a solid red line. 

(R/a) c ~ 5, corresponding to M » 360 particles, for ex- 
cess defects to proliferate in the ground state of a spher- 
ical crystal [22|. The precise value of (R/a) c depends on 
details of the microscopic potential but its origin is easily 
understood by considering just one of the 12 charge +1 
disclinations required by the topology of the sphere. In 
flat space such a topological defect has an associated en- 
ergy that grows quadratically with the size of the system, 
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since it is created by excising a 60° wedge of material and 
gluing the boundaries together |2^| . The elastic strain en- 
ergy associated with this defect grows as the area. In the 
case of the sphere the radius plays the role of the system 
size. As the radius increases, isolated disclinations be- 
come much more energetically costly. This elastic strain 
energy may be reduced by the formation of linear dislo- 
cation arrays, i.e. grain boundaries. The energy needed 
to create these additional dislocation arrays is propor- 
tional to a dislocation core energy E c and scales linearly 
with the system size. Such screening is inevitable in flat 
space (the plane) if one forces an extra disclination into 
the defect-free ground state. Unlike the situation in flat 
space, grain boundaries on the sphere can freely termi- 
nate [22, 13 EE |23, and our experimental results 
confirm these theoretical expectations. 

One systematic approach to determining the ground 
state of a collection of M particles distributed on the 
sphere and interacting via an arbitrary repulsive poten- 
tial [22I W\ treats the disclination defects themselves 
as the fundamental degrees of freedom, with the 6- 
coordinated particles forming a continuum elastic back- 
ground. The agreement between the predicted and ob- 
served values of (R/a) c supports the validity of this the- 
oretical approach. The original particle pair potential 
is replaced by a long range defect pair potential given 

by x{0) oc i? 2 (l + dz^), for a pair of defects 

separated by an angular distance j3. The potential is 
attractive for opposite charged defects and repulsive for 
like-charged defects. The underlying microscopic poten- 
tial enters only in determining the proportionality con- 
stant (equivalent to an elastic Young modulus) and E c . 
Many predictions of this model are therefore universal 
in the sense that they are insensitive to the exact micro- 
scopic potential. This enables us to make definite predic- 
tions even though the colloidal potential is not precisely 
known. It also means that our model system serves as 
a prototype for any analogous system with repulsive in- 
teractions and spherical geometry. To further test the 
validity of this approach, we show a typical ground state 
for large M in Fig. 03 The system size here is R/a = 12, 
similar to the droplet in Fig.QD. The results are remark- 
ably similar to the experimentally observed configuration 
in Fig.QjD; the only difference is a result of thermal fluc- 
tuations, which break the two defect scars in the exper- 
iment. This agreement between theory and experiment 
also provides convincing evidence that these scars are es- 
sential components of the equilibrium crystal structure 
on a sphere. 

The theory predicts that an isolated charge +1 discli- 
nation on a sphere is screened by a string of dislocations 
of length cos- 1 (5/6) J R » 0.59i? Q. We can use the vari- 
able linear density of dislocations to compute the total 
number of excess dislocations N in a scar. We find that 
N grows for large (R/a) as f [yTT - 5 cos -1 (5/6)] f « 




FIG. 3: Model grain boundaries. This image is obtained from 
a numerical minimization, based on the theory of 22], for a 
system size comparable to the droplet in Figs.^C,D). 

0.41(S) ; independently of the microscopic potential. 
This prediction is universal, and is in remarkable agree- 
ment with the experiment, as shown by the solid line in 
Fig. El 

We expect these scars to be widespread in nature. 
They should occur, and hence may be exploited, in 
sufficiently large viral protein capsids, giant spher- 
ical fuUerenes, spherical bacterial surface layers (s- 
layers) and the siliceous skeletons of spherical radiolaria 
(aulosphaera) [iij . provided that the spherical geometry 
is not too distorted. Terminating strings of heptagons 
and pentagons might serve as sites for chemical reactions 
or even initiation points for bacterial cell division |l2| and 
will surely influence the mechanical properties of spheri- 
cal crystalline shells. 
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